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In the course of its operation, the EGRET experiment detected high-energy 7-ray sources at energies above 
100 MeV over the whole sky. In this communication, we search for large-scale anisotropy patterns among 
the catalogued EGRET sources using an expansion in spherical harmonics, accounting for EGRET's highly 
non-uniform exposure. We find significant excess in the quadrupole and octopole moments. This is consistent 
with the hypothesis that, in addition to the galactic plane, a second mid-latitude (5° < |5| < 30°) population, 
perhaps associated with the Gould belt, contributes to the 7-ray flux above 100 MeV. 



1. Introduction 

Our rudimentary understanding of the GeV 7-ray sky was greatly advanced in 199 1 with the launch of the Ener- 
getic Gamma Ray Experiment Telescope (EGRET) on board the Compton Gamma Ray Observatory (CGRO). 
The science returns from EGRET observations have exceeded pre-launch expectations, increasing the number 
of known 7-ray sources from 1-2 dozen to the 271 listed in the third EGRET (3EG) catalog Q]. However, of 
this multitude of sources, only about half have been definitively associated with known astrophysical objects. 
One reason so few sources have been uniquely identified is due to a typical location uncertainty of up to 1° for 
an EGRET detected gamma-ray source, an area typically containing several potential candidate astrophysical 
objects. Consequently most of the 7-ray sky, as we currently understand it, consists of unidentified sources. 

The EGRET sources appear to arise from several distinct populations. There is an almost isotropic distribution 
of ~ 100 low/high variability objects, nearly all of which have now been plausibly associated with flaring 
blazars, a radio bright sub-class of active galactic nuclei |2||3|. The remaining sources are believed to be 
within the Milky Way, and there is some speculation that they may comprise two distinct classes. The first 
class consists of bright sources with generally hard spectra with galactic latitude |6| < 5° of the Galactic Plane, 
and are thought to be bright objects up to aproximately 10 kpc away. There is some evidence of a second 
population at mid-Galactic latitudes (5° < |6| < 30°) consisting of fainter objects with generally steeper 
spectra. It has been hypothesized that this population is associated with the Gould belt [4, 5J, an asymmetric 
structure in a great circle on the sky tilted 20° to the Galactic Plane. The Gould belt comprises massive O and 
B stars as well as clusters of molecular clouds and expanding interstellar gas at distance of around 200 pc from 
earth. The observation of spatial coincidence 1 6 , 7 8| of low-latitude (|5| < 10°) EGRET sources with massive 
OB stars and supernova remnants provides additional indication that the mid-latitude sources originate in the 
Gould belt. 

Figure[nshows the intensity map of EGRET sources. Though the correlation with the Galactic Plane is obvious, 
the hypthosized Gould belt population is not easy to extract using an eyeball fit. In this communication, we 
apply a power spectrum analysis (51 in an attempt to tease this component out of the bright Galactic Plane 
sources. 
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Figure 1. Intensity map of EGRET 7-ray sources with Figure 2. Two-dimensional EGRET detectability map for 
energies above 100 MeV. 7-i'ay sources (in units of cm^ s). 

2. EGRET Source Detectability 

For > 100 MeV, the significance for a detection of an isolated point source with EGRET scales as 

oc F^/S/B , (1) 

where F is the flux, £ is the exposure, and B is the intensity of the diffuse 7-ray emission at the region of 
the source flQ]. Following fl 1"|, we derive the detectability map considering these three ingredients in turn. 
First, the populations used in the Monte Carlo simulations are prepared such that they reproduce the observed 
distribution of the number of sources per flux bin, globally as well as locally. Second, it is necessary to account 
for the highly non-uniform exposure of the instrument. Thirdly, the diffuse 7-ray background is compensated 
using the EGRET diffuse model, corresponding to the commonly believed cosmic ray/gas distribution coupling 
(dynamic balance). 

With these considerations in mind, we construct the EGRET detectability map by determining 95% CL limits 
for a grid on the sky. We determine the limits using the so-called P1234 time span 1 1 1, i.e., we request that 
the maximally significant measurement of the EGRET source has been obtained in the P1234 joint viewing 
period. This tends to eliminate flaring sources from our sample. The map is shown in Figure|2] To obtain the 
detectability weights we normalize each angular bin to unity. 



3. Angular Power Spectrum 

The spherical harmonics form a useful set for expansions of the intensity over the celestial sphere. For a 
uniform sky coverage, 

00 t 

/(n)=^ ^ aemYem{n) , (2) 

1=0 m=-l 
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where n is a unit vector that denotes the direction in the sky of each source in the survey. ' To account for the 
nonuniform EGRET exposure, we define a weighted intensity 

1 ^ 1 

^(") = ^E-^(-'-.), (3) 

where tOj is the relative detectability at sky direction nj,J\f is the sum of the weights l^J^, and N is the number 
of sources. Since the eigenvalues of the Yf™ expansion are uniquely defined 



J I{n) y,„(n) d^n , (4) 



the replacement of Eq. (|3} into Eq. (0} leads to the explicit form of the coefficients for our set of = 124 
EGRET sources 

1 ^ 1 

a^m = 77 V — Yim{nj) . (5) 

yv — ; ojj 

J=l J 

The coordinate independent angular power spectrum of fluctuations, 

m = E aL , (6) 

provides a gross summary of the features present in the celestial distribution together with the characteristic 
angular scales. The power in mode £ is sensitive to variation over angular scales of radians. 

If the steady EGRET sources in fact originate from two overlapping great circles associated with the Galactic 
Plane and the Gould belt, one would expect a strong signal in the quadrupole (i = 2) moment describing the 
main component, and a fainter signal in the octopole {£ — 4) moment indicating the spread from a perfect 
quadrupole. Note that the Gould belt is tilted 20° relative to the Galactic Plane, so one would expect its 
presence to be mainifest in the octopole term which is sensitive to angular scales of about 15°. 

For £ > 1, the expectation for an isotropic distribution observed with uniform exposure is = {AnN)^^ , and 
the RMS fluctuations are (T2| 

AC, = V2/(2i?+l) Ci . (7) 

Figure |3l shows a comparison of the distributions of angular power spectra for isotropically distributed sources 
observed with uniform exposure (Ci) and with the EGRET exposure (Ci). It can be seen that for the four multi- 
pole moments of interest, the spread of the distributions corresponding to uniform and non-uniform exposures 
are similar 

Figure|4]compares the observed angular power spectrum of EGRET to the spectrum expected for an isotropic 
distribution. Clear excesses are seen for the quadrupole and octopole moments. The probability of the excess 
for either moment to be consistent with isotropy is estimated (from a fit to the distributions shown in Fig.|3} 
to be less than 0(10^^"). This supports the hypothesis that, in addition to the Galactic Plane, a mid-latitude 
population of sources, associated with the Gould belt, contributes to the observed 7-ray flux above 100 MeV. 



^Thi'oughout this work we use real- valued spherical hai'monics, which are obtained from the complex ones by substituting, e"" 

V2sin(m(/i),if m < 0, e*™"^ ^ V2cos(m0), if m > 0, ande""'^' 1 if m = 0. 



4 



Luis Anchordoqui et al. 




JO 



-2 

10 k 



+ random distribution 
® EGRET data 



3 4 

moment 



Figure 3. Comparison of the distributions of angular Figure 4. The angular power spectrum of the EGRET 
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